Cavity QED nondemolition measurement scheme using quantized atomic motion 
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Considering ultracold atoms traversing a high-Q Fabry-Perot cavity, we theoretically demonstrate 
a quantum nondemolition measurement of the photon number. This fully quantum mechanical 
approach may be understood utilizing concepts as effective mass and group velocity of the atom. 
The various photon numbers induce a splitting of the atomic wave packet, and a time-of-flight 
measurement of the atom thereby reveals the photon number. While repeated atomic measurements 
increase the efhciency of the protocol, it is shown that by considering long interaction times only a 
few atoms are needed to resolve the photon number with almost perfect accuracy. 
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I. INTRODUCTION 

Any measurement of a quantity A of a quantum sys- 
tem has an impact on the system itself. For example, 
exactly determining A implies that any knowledge of a 
conjugate variable is lost. Moreover, the detection might 
even be destructive by nature, e.g. the standard way of 
measuring the electromagnetic field is by photocounting 
detectors where photons are actually absorbed. However, 
there exist situations where a measurement of some quan- 
tity A does not induce any back-action quantum noise on 
A. That is, provided A is a constant of motion, then a 
second measurement of A, after some time-delay relative 
the first measurement, would reveal the same value as 
obtained in the first measurement. This is called a quan- 
tum nondemolition (QND) measurement 

By now, QND measurements of the optical field have 
been demonstratedjsoth for an optical fiber 0] and in 
cavity QED setups [3j]. In the cavity experiments to date, 
an atomic Ramsey interferometer technique has been uti- 
lized. The different phases aquired for an atom interact- 
ing dispersively with the cavity field either in its excited 
state |e) or ground state \g) contain information about 
the photon number n. Repeating the atomic QND mea- 
surements sufficiently many times determines the photon 
number definitely. 

The kinetic energy of the atoms in these experiment 
greatly exceeds the atom-field interaction energy, and 
thereby they can be safely treated by classical means. For 
ultracold atoms on the other hand, mechanical effects in- 
duced by the light fields become important. Such actions 
are indeed the building blocks for cooling and trapping 
of neutral atoms In the cavity QED community, it 
has long been known that the atom-field dynamics may 
be considerably modified by treating the atomic motion 
quantum mechanically together with taking spatial mode 
variations into account [5| . The system now includes ad- 
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ditional degrees of freedom that become correlated with 
the cavity field. Already back in 1989 it was shown that 
deflection of a beam of ultracold atoms resonantly inter- 
acting and transversely scattered from a quantized stand- 
ing wave field will depend on the actual photon distribu- 
tion of the field @. For a classical field, the same effect, 
named optical Stern-Gerlach, was presented in Ref. [7|. 
These observations led to the idea of performing QND 
measurments of the cavity field using ultracold atoms 
dispersively interacting with the field. Reference |8|| con- 
sidered a beam of ultracold atoms transversely passing 
a Fabry-Perot cavity. In the Raman-Nath regime, valid 
for very short interaction times, it was shown that the 
deflection of the atomic beam after interacting with the 
cavity field depends on the number of photons. Thus, 
by recording the atomic positions of a sequence of atoms 
having interacted with the cavity field, a QND measure- 
ment of the photon number is possible. 



In this paper we address a different QND measure- 
ment where the time-of-flight of the atom is recorded. 
Instead of studying the deflection of transversely pass- 
ing atoms, we consider atoms traversing the Fabry-Perot 
cavity along its axial direction. Our treatment is fully 
quantum mechanical, taking into account for; the atomic 
scattering effects occurring as the atom enters and exits 
the cavity, the quantum pressure arising from the atomic 
kinetic energy term (hence going beyond the Raman- 
Nath regime) , as well as the uncertainty of the velocity 
of the incoming atom. The scheme is first studied in a 
semi- analytical model which relies on the concepts of ef- 
fect mass and group velocity for the atom. By imposing 
a sort of single-band approximation, the semi-analytical 
model demonstrates that a single atom is sufficient for 
performing the QND measurement in the ideal situation 
provided the effective interaction time is long enough. 
The full system, going beyond single-band and taking 
into account for a finite cavity, is considered numerically, 
and it is found that typically around ten atoms is enough 
for achieving a highly efficient QND measurement. 
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II. MODEL SYSTEM 



III. SEMI-ANALYTICAL ANALYSIS 



We consider an ultracold two-level atom sent through 
a Fabry-Perot cavity along its axial axis. The internal 
ground and excited atomic states are labeled \g) and |e), 
and their energy difference hTl. The atom interacts with 
a single cavity mode with frequency uj. Moving to an in- 
teraction frame and imposing the rotating wave approx- 
imation, the Hamiltonian reads 



H' 



2m 



hA 



g{X)a^a- +g*iX)a+a . (1) 



Here, X and P are the atomic center-of-mass position 
and momentum respectively, A = il — w is the atom-field 
detuning, g(X) the effective position dependent coupling, 
(d) the creation (annihilation) operators for the field, 
and the Pauli-operators are az = |e)(e| — \g)(g\, (t^ — 
\e){g\, and a~ — \g){e\. For a cavity of length L we have 



In the previous section we introduced the model Hamil- 
tonian. Furthermore, it was assumed that L 3> 27rfc~^, 
making sure that the system is quasi periodic. Neglect- 
ing boundary effects arising from having a finite cavity 
length, the corresponding eigenvalue problem relaxes to 
the Mathiue equation 



(6) 



As a periodic operator, the spectrum E"{k) for given 
n is characterized by a band index v = 1,2,3, ... and a 
quasi momentum q defined within the first Brillouin zone 
(—1 < g < 1). The corresponding eigenstates 
are the so called Bloch functions. A typical spectrum 
is envisaged in Fig. [1] The figure shows the first three 
energy Bloch bands 
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The length L of the cavity is assumed to be much larger 
than the wavelength of the field, and hence, the Hamil- 
tonian is quasi-periodic, k is the wave number and go 
the effective atom-field coupling. Letting £^2r = de- 
fine the characteristic energy, we scale the parameters 
accordingly 



go = go/E2r, A = A/i?2r, 
t = tE2r/fi, X — kX, 



(3) 



where t is the unsealed time. Using the fact that within 
the rotating wave approximation, the number of exci- 
tations N — a^a + (Jz/^. is preserved, the Hamiltonian 
may be written on block form within the states \n^g) 
and |n — l,e), \n) being the Fock state with n pho- 
tons. Thus, we have H' = Hq (E) H'l (E) H2 ® ■■■ where 
-dV2dx2 _ A/2 and 



H'o 



H' 



2rfi2 



A 

2 



go cos{x)^/n 



go cos(x)y/n 



A 
' 2 



(4) 



Here we have taken go to be real. 

For ultracold atoms in the dispersive regime, A 3> 
go\/n, we adiabatically eliminate the excited atomic level 
|e) to obtain a Hamiltonian describing the dynamics for 
the atomic ground state alone. Following standard pro- 
cedures 0, one derives 



2dx^ 



(5) 



where C/ = ^ is the amplitude of the single photon dipole 
induced potential. 




FIG. 1: The first three Bloch bands E^{q) {v = 1, 2, 3). The 
amplitude Un — 0.5. AH parameters are dimensionless. 



Typically, incoming atoms are not in pure momentum 
eigenstates. The atomic velocity selection is never perfect 
causing an uncertainty in the mean momentum. This is 
taken into account by considering initial Gaussian states 
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where Ap is its width determined by the velocity uncer- 
tainty in the state preparation, and po the mean initial 
momentum. The time evolution of this state is rendered 
by the Hamiltonian ([5]). For a small coupling Un and 
moderate spreading Ap <ti 1, the atomic state will pre- 
dominantly populate a single Bloch band v' with average 
quasi momentum go = po — v' . By expanding the corre- 
sponding energy around qo 



E:,{q)^E{qo) + v'l{qo){q-qo) + T; 
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where = dE'^,{q)/dq\q=q„ and to,; (go) = 

{(fEl^,{q)/dq^y it follows that, for a given n, the 
time evolved state in ^-representation approximates [lol | 



where 
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and Vg = Vg{qQ) and to^ = TOj^(go)- Within these ap- 
proximations, the wave packet preserves its Gaussian 
form, moves with a group velocity Vg, and spreads ac- 
cording to the effective mass to^. 

For a general initial state of the cavity field |(/)) — 
^^QC„|n), we obtain the single-band approximated 
time evolved atom-field state 
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FIG. 2: The atomic density HI} at time t = 400. The field is 
initially in a coherent state with n = 4 and the inset numbers 
indicate the corresponding number of photons. The param- 
eters are po = 2.58, Ap = po/50, and U — 0.7. The initial 
velocity po implies that the atomic wave packet evolves on the 
third Bloch band. 



and the corresponding atomic density 



Pat{x,t) 
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Due to the uncertainty of photon numbers, it follows that 
the atomic state will split into a set of Gaussians. Hence, 
whenever the distances between consecutive Gaussian 
wave packets exceeds their widths, a measurement of the 
atomic position will reveal the photon number n. Equiv- 
alent to a position measurement is a time-of-flight mea- 
surement, where the time it takes the atom to traverse 
the cavity is recorded. The atomic density Pat{x, t — 400) 
is depicted in Fig. O The initial cavity field is a coherent 
state 



l"lV2. 
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with an average number of photons n = jap = 4. For 
the current set of parameters, the wave packet evolves 
on the third Bloch band. The inset numbers give the 
corresponding number of photons n. For n = 4, the 
population of the Fock state with n = 10 is less than one 
percent and therefore only the first ten Fock states are 
seen. 

For the example of Fig. [51 a single atomic detection 
is most likely sufficient for determining the photon num- 
ber. However, for shorter interaction times the Gaus- 
sians overlap and a single measuremnt cannot resolve 
the photon number. Nonetheless, by repeating the mea- 
surement procedure for a second and third atom and so 
on, the photon distribution will finally pick a single Fock 
state \n) despite the non-zero Gaussian overlaps. Nat- 
urally, the less resolved the peaks are the more atoms 



are needed. Stated in other words, a strong atom-field 
correlation implies fewer atomic detections. Consider- 
ing pure initial states and a closed system, we employ 
the von Neumann entropy as an estimate of the amount 
of atom-field correlation/entanglement. For the reduced 
field density operator Pf{t) = Trat[p(x, i)], where p{x,t) 
is the full system density operator and the trace is over 
the atomic motional degrees of freedom, we have the von 
Neumann entropy 



S}{t)^-TTf[p^{t)Hpi[t)]]. 



(14) 



The trace is over field degrees of freedom. The maximum 
entropy is given by S'j""'^) = - S„ l^nP ln[|c„p], with c„ 
the initial photon amplitudes. Using the same parame- 
ters as in Fig. [21 we present the corresponding entropy in 
Fig. [31 The dashed line gives the maximum entropy, and 
it is clear that Sf{t) approaches this value in the long 
time limit. On the other hand, for times t less than say 
100, a series of atomic detection is most likely required 
for an efficient QND measurement. 



IV. NUMERICAL ANALYSIS 

A. Complications due to a realistic system 

In realistic experimental setups, several complicating 
effects arise that were not fully addressed in the previous 
section. In this section we will take them into account 
numerically. 

By entering the cavity, we assume that the atom feels 
a fairly sudden turn on of the dipole induced cavity po- 
tential. To model the potential Q for the whole x-axis, 
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FIG. 3: Time evolution of the von Neumann entropy (|14p . 
Parameters are the same as in Fig. (2] The amount of entan- 
glement is rapidly increasing and asymptotically reaches its 
maximum value 5* 
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indicated by the dashed line. 



we introduce an envelope function such that 



, , x-Ll2\ , fx+Ll2 
tanh — — — — tanli ' 
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where Xs determines the slope of the envelope function 
around x = ±L/2 and we naturally chose <C L. The 
rapid change in V{x) in the vicinity of a; = ±L/2 will 
induce some backward scattering of both the incoming 
and outgoing atomic wave packets. Thereby, especially 
for large amplitudes Un, part of the wave packet will not 
reach the detector. Experimentally, backward scattered 
atoms may be ignored and we will focus only on the for- 
ward scattered atoms, those who are detected. The mea- 
sured time tn is presumably the time between atomic 
state preparation and the detection. As no backwarded 
scattered atoms are recorded, t„ includes the free space 
time propagation tfr plus the time t„ spent inside the 
cavity; tn = Tn + tfr- Since the atom- field interaction is 
dispersive, the atomic velocity before and after the cavity 
are the same and independent of n, and consequently t„ is 
an indirect measure of the group velocity Vg{qo) « L/tu- 
Yet another complication is the fact that the true state 
of the atom ([7]) is not restricted to a single Bloch band 
I/'. From the symmetry of the Hamiltonian, and given 
n, it follows that the atomic state expressed in Bloch 
functions is 
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d:i;{q,0)\rAQ))dq. (16) 



Here, ip{q, 0) is ip{p, 0) with the real momentum replaced 
by the quasi momentum [ll[, and the d" are the proper 
weights given Un. For small Un, the coefficient w 1 
with v' = pq — qo. For the parameters employed in 
the previous section one finds that all jd", | is noticeable 



smaller than 1, and the atomic wave packet will therefore 
split up into sub-packets as it enters the cavity. The dif- 
ferent parts deriving from atomic propagation according 
to the corresponding Bloch bands. 

It is clear that it is a trade-off between having a large 
separation of group velocities , and small atomic scat- 
tering and splitting. The first favors large amplitudes 
[/, while effects originating from scattering and splitting 
are decreased for small values of U. Despite this, we will 
now demonstrate that by repeated measurements, the ef- 
ficiency can be made asymptotically close to unity. 



B. Results 

The Schrodinger equation ^ is solved using the split- 
operator method |12l |. The atmic intial condition is a 
Gaussian with width A^; = 15 S> 27r/c~^, initial position 
xq = — L/2 — 70, and initial momentum pq = 3.75. The 
cavity length L is varied, while Xg — 0.2 and U = 0.7 
are kept fixed throughout. The initial cavity field is as 
before, i.e. a coherent state with n — 4. Note that 
the atomic wave packet is initialized in the regime where 
the amplitude of the cavity field is approximately zero; 
V{—L/2 — 70) w 0. The time propagation tf is per- 
formed till the various fowardly scattered atomic wave 
packet components have left the interaction region. At 
this instant t/, due to the back scattering we renormalize 
Pat{x,tf) accordingly 
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Pat(x,tf) 



X < L/2 
X > L/2. 



(17) 



In Fig. [¥] (a) we present the renormalized atomic den- 
sity for X > L/2. The cavity length for this example is 
fairly long, L = 1400. Compared with the ideal situation 
of Fig. [21 we note that in this more realistic situation the 
various atomic wave packet components ipn{x,tf) are less 
resolved from each other. In Fig.[3] (b) we give the corre- 
sponding components ^^(x,t/) renormalized within the 
interval x £ [L/2, oo]. The numbers represent the photon 
number n. The peaks around x w 2000, corresponding 
to photon numbers n = 6, 7, 8, originate from the wave 
packet component propagating on higher excited Bloch 
bands. Expectedly, this component becomes significantly 
populated only for strong field amplitudes Un. 

For the example of Fig. [31 a single atomic measure- 
ment will not reveal the photon number with very high 
efficiency. However, as mentioned in the previous section, 
repeated atomic measurements will improve the scheme 
considerably. First we note that numerically, instead 
of considering a time measurement we can equally well 
freeze the time evolution and make a position measure- 
ment. After the first atom has been recorded with a cor- 
responding position Xr, the photon distribution becomes 



|c«|2 = |v;(a;„i/)|2|c„|ViV, where TV = EJcJ.^ 
the second atom traversing the cavity, its atomic den- 
sity p}p{x,t) is given by Eq. (jT^j) with c„ replaced by 
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FIG. 4: The renormalized atomic density Pat{x,tf) (a) and 
the different atomic wave packet components tp'„{x,tj) (b). 
The cavity length L = 1400 and final time tf = 660. Other 
dimensionless parameters 15, xo = -770, po = 3.75, 

n = A,U = 0.7, and Xs = 0.2. 



c„ . Each atomic measurement acts as a "photon fil- 
ter" [isl ]. modifying the photon distribution with the 
weights tp'^{xr,tf ). 

For the numerical simulation, the position Xr of atom j 
is randomly picked according to the corresponding prob- 
ability distribution p'^p{x,tf). Once the position has 
been determined, the photon distribution is adjusted ac- 
cordingly and a new atomic density for atom j + 1 is 
calculated. The process is repeated until only a single 
photon component survives the filtering. The number of 
iterations needed depends on; the randomly picked num- 
bers and on how well separated the atomic wave packet 
components are. In general, large L implies less atomic 
measurements . 

The results of three simulations for different cavity 
lengths L are presented in Fig. [5l In the first case (a), 
the parameters are the same as for Fig. |31 Already after 
three atomic measurements the photon distribution has 
collapsed to approximately a single Fock state. For the 
cases with L = 600 (b) and L — 200 (c), considerably 
more atoms are needed in order to single out a photon 
number. 

Each atomic measurement is projective, leaving the 
field in a pure state = J2n '^'n\xr,tf)cn'' \n) /Nj, Nj 
being the proper normalization constant. The filtering 
projection onto a single Fock state seen in Fig. [5] can also 
be demonstrated via the phase space distributions, e.g. 
the Husimi Q-function [ij] 

QW(«) = i*i^. (18) 

TT 

Here, |a) is a coherent state with amplitude a. For a 
coherent state, the Q-function is Gaussian, while for a 
Fock state it is a circle with radius |a| = ^/n. Figure [5] 




FIG. 5: Simulations of the QND measurement using succes- 
sive atomic measurements. The three plots correspond to 
L = 1400, L = 600, and L = 200 respectively. For longer 
cavities, less atomic measurements are required to have an 
efficient QND measurement. The rest of the parameters are 
the same as in Fig. |4]except xo ~ — L/2 — 70. 



gives four examples of for j = (a), j = 1 (b), 

j — 5 (c), and j = 10 (d). The initial state and the 
parameters are the same as those of Fig. [5] (b). In plots 
(b) and (c), there are two photon numbers dominating, 
n = 4, 5. In (d), almost all population resides in the 
n = 3 state. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

The idea behind our scheme is different from that of 
Ref. Q. In due to virtual exchange of photons with 
the cavity field, atoms are deflected perpendicularly with 
respect to their initial velocity. A position measurement 
is therefore an indirect measure of the number of pho- 
tons that has been exchanged. Even for short interac- 
tion times (imposing the Raman-Nath approximation), 
the atom can acquire a certain number of 2Hk momen- 
tum kicks, where the factor 2 comes from the fact that 
the interaction is dispersive and absorption of a photon 
by the atom is always accompanied by emittance of a 
photon. The measurement then projects onto any of the 
momentum eigenstates p = 2rHk, r = 0, 1, 2, .... In 
the present scheme, on the other hand, the time-of-flight 
measurement indirectly gives the wave packet velocity 
which for weak couplings approximate the group velocity 
Vg , characterizing the average velocity (p) / m for the cor- 
responding Bloch state. To make these arguments more 
transparent, in Fig. [7] we display the time evolution of 
(p) for the L = 600 cavity and n = 0, 1, 2, 3, 4. The 
plot makes clear that the velocity decreases for increas- 
ing photon numbers, and also that the atom regains its 
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FIG. 6: Evolution of the Q-function (|18p after the positions 
of j = 0, 1, 5, 10 atoms have been recorded. For the ini- 
tial state (a), the Q-function is a Gaussian centered around 
{ar,at) = (2,0). Even after 1 atomic measurement (b), the 
phase space quasi distribution does show a circular structure 
characterizing the Fock number state. After 5 measurements 
(c), the circular structure has actually declined, while after 10 
atoms (d) the distribution is almost perfectly circular with a 
radius |q:| « 1.7 corresponding to the n — 3 Fock state. The 
parameters are as in Fig. [S] (b), but the sequence of random 
numbers Xr is not the same as those used for that plot. 



initial velocity after exiting the interaction region. 




The measurement efficiency is enhanced by consider- 
ing long cavities. However, long cavities naturally im- 
plies long interaction times and cavity losses may be- 
come significant. For short cavities, the individual atomic 
interaction times are shorter but, on the other hand, 
more atomic measurements are needed. This again might 
cause long total operational times. One way of de- 
creasing the total process time is by using feedback- 
techniques [IBl • Nonetheless, even when cavity losses be- 
come important the scheme can be useful. For a lossy 
cavity, one typically considers an external pumping of 
the cavity, keeping the field in a coherent state with am- 
plitudes determined from balancing the pump and loss 
rates. If the time-scale for the field is much shorter than 
that of repeated atomic measurements, the field attains 
its steady state between each measurement. The result 
of the measurements will then reveal the steady state 
photon distribution. 

As a summary, we have introduced a QND measure- 
ment scheme of the photon numbers in a cavity. It 
relies on time-of-flight measurements of atoms travers- 
ing a Fabry-Perot cavity along its axial direction. The 
field intensity n directly affects the atomic velocity while 
traversing the cavity; the velocity drop is increased for 
larger photon numbers n. This can be explained us- 
ing the language of group velocities for particles mov- 
ing in periodic potentials. Employing this concept, we 
argued that for a sufficiently long cavity a single atom 
can resolve the photon number non-destructively. By 
decreasing the cavity length, repeated atomic measure- 
ments are required. Our analytical results were numer- 
ically verified considering more realistic cavity QED se- 
tups. In particular, we found that typically more than 
single atomic measurements are needed for reliable QND 
measurements. We further explained that even for lossy 
cavities the scheme is of interest as it can be used as 
a tool to measure the full steady state photon distribu- 
tion. Getting access to the full state of the field would 
require more sophisticated approaches like tomography 
measurements 
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